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ABSTRACT 

Supernova (SN) explosions enrich the intracluster medium (ICM) both by creating and dispersing metals. We 
introduce a method to measure the number of SNe and relative contribution of Type la supernovae (SNe la) and 
core-collapse supernovae (SNe cc) by directly fitting X-ray spectral observations. The method has been implemented 
as an XSPEC model called snapec. snapec utilizes a single-temperature thermal plasma code ( apec ) to model the 
spectral emission based on metal abundances calculated using the latest SN yields from SN la and SN cc explosion 
models. This approach provides a self-consistent single set of uncertainties on the total number of SN explosions 
and relative fraction of SN types in the ICM over the cluster lifetime by directly allowing these parameters to 
be determined by SN yields provided by simulations. We apply our approach to XMM-Newton European Photon 
Imaging Camera (EPIC), Reflection Grating Spectrometer (RGS), and 200 ks simulated Astro-H observations of a 
cooling flow cluster, A31 12. We find that various sets of SN yields present in the literature produce an acceptable fit 
to the EPIC and RGS spectra of A3112. We infer that 30.3% ± 5.4% to 37.1% ± 7.1% of the total SN explosions 
are SNe la, and the total number of SN explosions required to create the observed metals is in the range of 
(1.06 ± 0.34) x 10 9 to (1.28 ± 0.43) x 10 9 , from snapec fits to RGS spectra. These values may be compared to the 
enrichment expected based on well-established empirically measured SN rates per star formed. The proportions 
of SNe la and SNe cc inferred to have enriched the ICM in the inner 52 kpc of A31 12 is consistent with these 
specific rates, if one applies a correction for the metals locked up in stars. At the same time, the inferred level of 
SN enrichment corresponds to a star-to-gas mass ratio that is several times greater than the 10% estimated globally 
for clusters in the A31 12 mass range. 
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1. INTRODUCTION 

The use of Type la supernovae (SNe la) as standardizable 

candles in the discovery of dark energy (Riess et al. 1998; 
Perlmutter et al. 1999) installed the nature of SN la progen- 
itors and the physics of the explosion as one of the prime 
problems in astrophysics. The evolution of the SN la rate is 

one of the promising methods for unveiling SN la progenitors 

(Ruiz-Lapuente & Canal 1998; Yungelson & Livio 2000). Mea- 

suring the evolution of the supernova (SN) rate is particularly 
important at higher redshifts where direct SN rate constraints 
are limited (Gal- Yam et al. 2002; Barbary et al. 2012). 

Galaxy clusters represent the largest scales of organized 
matter in the universe, making them unique laboratories for 
chemical enrichment of the universe from all possible sources 
including SNe la. The intracluster medium (ICM) emits X-rays 
due to the highly ionized gas which has been heated by infall 
from the intergalactic medium (Gunn & Gott 1972). The metals 
produced, during all stages of cluster formation and evolution, 
by stars and galaxies enter the ICM via SN explosions (De 
Young 1978) and strong galactic winds (Mathews & Baker 

1971). They remain in the ICM due to the cluster’s deep potential 
well. Since X-ray spectroscopy yields accurate measurements 
of metal abundances in the ICM, the large reservoir of metals 
in clusters of galaxies provides a unique way to probe SN rates 
and therefore SN la progenitor models on a universal scale. 

The launch of ASCA triggered studies of the relative con- 
tributions of SNe la and core-collapse supernovae (SNe cc) to 


ICM metal enrichment. These studies used specific elemental 
abundance ratios (e.g., nickel to iron) measured in clusters to 
distinguish between SN la models (Baumgartner et al. 2005; 
Dupke & White 2000; Dupke & Arnaud 2001; Mushotzky & 
Loewenstein 1997). With the launch of observatories with better 
spectral and spatial resolution, such as XMM-Newton, Chandra, 
and Suzaku, radial profiles of specific elemental abundances 
were used to determine the relative contribution of SNe la to the 
enrichment of clusters as a function of position (e.g., Matsushita 
et al. 2007; Sato et al. 2007; Million et al. 2011). High-resolution 
XMM-Newton Reflection Grating Spectrometer (RGS) measure- 
ments of specific elemental abundance ratios within the ICM 
have also been used to determine the integral yield of all the SN 
la and SN cc explosions and to distinguish between SN models 
with different levels of pre-enrichment of the progenitors and 
with different initial mass functions (IMF; Werner et al. 2006b; 
de Plaa et al. 2007; Grange et al. 2011). 

In this paper we introduce a new approach to determine the 
total number of SN explosions, and the relative contributions 
of different SN types, to the total enrichment of the ICM 
using high-resolution X-ray observations. The method relies 
on using the full set of nucleosynthetic yields produced by 
recent SN models for elements with emission lines in the 
X-ray bandpass to probe the total number of SN explosions 
and relative contribution of SNe la. We validate this approach 
using CCD resolution XMM-Newton European Photon Imaging 
Camera (EPIC) and high-resolution RGS observations of the 
cooling flow cluster A31 12. A brief summary of the method is 
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given in Section 2. The application of the method to analysis of 
XMM-Newton RGS, EPIC, and simulated Astro-H Soft X-Ray 
Spectrometer (SXS) data is outlined in Section 3. We provide 
our discussions and conclusions in Sections 4 and 5. 

2. XSPEC SUPERNOVA ABUNDANCE MODEL ( snapec ) 

2.1. A Brief Description of snapec 

Most of the metals from oxygen through nickel residing in the 
ICM were synthesized by SN la or SN cc explosions. We provide 
a new XSPEC model called snapec to determine the total number 
of SN explosions and relative contributions of SN types to the 
metal content of the ICM using elemental yields provided by 
SN nucleosynthesis calculations in the literature. The expected 
elemental yields of SNe la and SNe cc from massive ( 10-50 M 0 ) 
progenitor stars of various metallicity (0-1 A 0 ) have been 
intensively studied by several authors (Woosley & Weaver 1995; 
Iwamoto et al. 1999; Nomoto et al. 2006). For instance, SN 
models predict that SNe la produce significant amounts of iron 
(Fe), nickel (Ni), and silicon (Si), while, for SNe cc, large 
quantities of oxygen (O), neon (Ne), and magnesium (Mg) are 
produced but very little Fe and Ni escapes the compact remnant. 

snapec calculates the mass of the /th element (M? Ne ) in terms 
of the number of SN la (/V Ia ) and SN cc (N cc ) explosions that 
enrich the ICM and the yields per SNe la (yf) and SNe cc (y? c ), 


M SNe = yf + N cc (yfcj 

(1) 

= /V SNe (l + R)-‘ [Ry, Ia + [yf)]. 

(2) 


where R is the ratio of SNe la to SNe cc (R = /V la //V cc ), and 
/V SNc is the total number of SN explosions /V SNc = /V la + N cc . 
The SN yields yf and yf are obtained from published SN mod- 
els and stored in the snapec database in solar units. The database 
includes yields of 30 elements from SN la explosions of slow 
deflagration (W7, C-DEF), delayed detonation (WDD, CDDT, 
ODDT), and core degenerate scenarios (CDD) (Tsujimoto et al. 
1995, hereafter T95; Iwamoto et al. 1999, hereafter 199; Maeda 
et al. 2010, hereafter M10), and from SN cc yields from an 
extensive range of progenitor masses ( 10-50 M 0 ) and metal- 
licities (0-1 times solar) (Tsujimoto et al. 1995; Iwamoto et al. 
1999; Woosley & Weaver 1995, hereafter WW95; Nomoto et al. 
2006, hereafter N06). snapec adopts solar abundances for met- 
als with atomic number fl (e.g., helium, lithium, beryllium, 
boron, carbon, and nitrogen). 

The snapec model allows a selection of solar abundance 
standards — the database incorporates the sets of Anders & 
Grevesse(1989), Lodders (2003), or Asplundetal. (2009). Since 
the majority of hydrogen and helium in the ICM was produced 
during the big bang, the standard hot big bang nucleosynthesis 
(BBN) model predicts the primordial helium abundance by unit 
mass Yp to be 0.2565 ± 0.0010 (Izotov & Thuan 2010) which 
corresponds to 8.3% helium atoms in number assuming 0.3 solar 
abundances for metals (e.g., Markevitch 2007). The Asplund 
et al. (2009) solar abundances we use here predict 8.5% of 
helium which is consistent with BBN predictions. 

The SN cc yields of the /th element, averaged over the IMF 
(0), is 

/ cc \_ 

[y ‘ ' C<Km)dm ’ 

where m\ and mi are the lower and upper limits for the masses 
of SN cc progenitors. Here we assume that stars with masses 


M ^ 10 Mq and M A 50 37 do not explode as SNe cc. In this 
work we use the standard Salpeter IMF slope above 10 Mq, 
f = m ~ 2 35 (Salpeter 1955). 

It is straightforward to transform Equation (2) to an expres- 
sion for the abundance, by number relative to hydrogen, of the 
/th element. As long as the total mass fraction of metals is small 
and the He abundance is solar, 


Zf Mf 

(4) 

N SNe 

0 [ R yf + (y, cc )l . 

Mf( 1 + R) 1 ' ' ‘ /J 

(5) 


Mf is the mass of the /th element for an ICM with solar 
abundances, i.e.. 


Mf 


MiCMZf A i 


m N 


(6) 


where Zf is the abundance of the /th element in solar units, A, 
is the atomic mass of the /th element, and m,y = 2.27 x 10~ 24 
g/H-atom is the average nucleon mass per hydrogen atom 
(Asplund et al. 2009). 

This XSPEC model, snapec, has five free parameters: gas 
temperature (kT e ), the total number of SN explosions (/V SNe ) 
and the SN ratio (R) that enriches the ICM in the spectral 
region being analyzed, the cluster’s redshift, and the XSPEC 
normalization factor identical to those defined for other thermal 
plasma models (see below). In addition, there are two fixed 
parameters, numerical indices that specify the SN la and SN 
cc models from which the yields are drawn. The predicted 
spectrum is constructed using the single-temperature thermal 
model ( apec ) (Smith et al. 2001) and AtomDB 2.0.1 (A. 
Foster et al. 2012, in preparation) to determine the continuum 
and line emission for the metal abundances calculated using 
Equation (5). This method has the advantage of simultaneously 
providing a single set of uncertainties on the SN ratio and the 
total number of SN explosions that enrich the ICM since the 
formation of the cluster. The results presented in this work are 
obtained using the Asplund et al. (2009) solar abundances. 

snapec redefines the quantities in Equation (5) to express 
the abundances in terms the total number of SNe per 10 12 M 0 
of ICM plasma — i.e., rescaled to yield values appropriate for 
cluster cores. The total number of SNe is derived by multiplying 
by the ICM mass (in units of 10 12 M 0 ) within the spectral 
extraction region. This may be estimated from the XSPEC 
normalization, which is proportional to the emission measure. 


N = 


10 ~ 14 

4tt (1 + z) 2 D\ 


J n e tindV, 


(7) 


where n e is the electron number density, l) A is the angular 
diameter distance, and z is the cluster’s redshift. 


3. APPLICATIONS TO X-RAY SPECTRA 

Accurately determining the elemental abundances is crucial 
for understanding the relative contributions of different SN types 
to the metal enrichment of galaxies and intracluster gas. Earlier 
attempts to determine the relative fraction of SNe la using 
XMM-Newton measurements of individual elements such as Si, 
S, Ar, Ca, Fe, and Ni in the ICM (Werner et al. 2006b; de 


2 


The Astrophysical Journal, 753:54 (9pp), 2012 July 1 


Bulbul, Smith, & Loewenstein 


Table 1 

XMM-Newton RGS Observations of A31 12 


Obs. ID 

R.A. 

Decl. 

Exposure 

Clean 




Time 

Time 


(J2000) 

(J2000) 

(ks) 

(ks) 

0603050101 

03 17 57.4 

-44 14 14.7 

118.9 

100.1 

0603050201 

03 17 57.4 

-44 14 12.8 

80.8 

80.4 


Plaa et al. 2007) concluded that the T95 and 199 SN models 
overestimated Ar and Ca abundances, and could not produce 
acceptable fits to XMM-Newton data, de Plaa et al. (2007) also 
inferred a number ratio of SNe cc to SNe la (using W7 yields) 
of ~3.5, based on XMM-Newton observations of 22 clusters. 
Grange et al. (2011) studied EPIC and RGS observations of 
two groups of galaxies, NGC 5044 and NGC 5813, and found 
that ~30%-40% of SNe contributing to metal enrichment of 
the intragalactic gas were Type la while the other ~60%-70% 
of SNe then were core collapse. This result was consistent 
with what was previously found for M87 by Werner et al. 
(2006a). Simionescu et al. (2009) reported that WDD3 or W7 
SN la models are needed in order to reproduce the low Si 
abundance in Hydra A clusters of galaxies. They also noted 
that abundance patterns in the XMM-Newton observation of all 
other clusters preferred WDD1 and WDD2 models. However, 
these studies proceeded by first determining the best-fit values 
and uncertainties on individual elemental abundances, followed 
separately by fitting abundance ratios to determine the relative 
fractions of SNe la and SNe cc using the T95 and 199 models. 
By contrast, in this work we apply the yields obtained from 
SN models to immediately predict the spectrum, skipping the 
intermediate step required in previous studies. The best-fit R 
and /V SNc parameters, and parameter confidence levels, are 
thus directly estimated from X-ray spectra and their statistical 
uncertainties for a given set of SN yields, allowing one to use 
the goodness of fit to assess the accuracy of the predictions and 
compare different yield sets on a well-defined statistical basis. 


Subsystem (SAS) software. For details about RGS data reduc- 
tion and analysis, see Bulbul et al. (2012, hereafter Paper I). We 
fit the combined RGS first- and second-order spectra using an 
absorbed snapec model in the 7-22 A and 7-16 A wavelength 
range, respectively. The absorption column density was fixed 
to the Leiden/ Argentine/Bonn Galactic value (Kalberla et al. 
2005). The /V SNc , R. redshift, kT e , and normalization parame- 
ters were allowed to freely vary. The best-fit model parameters 
obtained from the snapec fits to RGS spectra for various SN 
models are shown in Table 2. Since, as in Paper I, fitting is con- 
ducted using the C-statistic and the C-statistic does not provide 
a direct estimate of the goodness of fit, we asses the goodness of 
fit using the corresponding y 1 value. These final y 2 values, ob- 
tained from the best fit determined by C-statistics, are reported 
in Table 2. 

We calculated the projected gas mass of A3 1 12 confined in the 
central 38" region corresponding to the effective RGS aperture 
(Paper I) based on the best-fit normalization parameter in the 
single-temperature thermal model ( apec ) fit to the RGS spectra 
using Equation (7). The apec normalization of 1 .28 x 1 0 2 gives 
a projected gas mass of (4.80 ± 0.70) x 10 11 M e within the 
38" (~52 kpc) region. We also calculated a spherical gas mass 
of (4.06 ± 0.43) x 10 11 M q within the 38" (~52 kpc) which 
is consistent with the projected mass obtained from the apec 
normalization. This spherical gas mass was calculated using the 
deprojection analysis and modeling of the XMM-Newton EPIC 
observations of A3112 (see Bulbul et al. 2010, 2012). Thus, 
a conversion factor 0.41 is applied to the SNe per 10 12 M G 
derived from spectral fits to obtain the total number of SNe, 
/V SNe parameters, which are also reported in Table 2. 

We also convert the R parameter obtained from snapec fits 
into the fractional contribution of SNe la to the total number of 
SN explosions as follows: 


(SNe la + SNe cc) 


3.1. XMM-Newton RGS Obser\>ations of A3 112 

XMM-Newton' s high spectral resolution RGS instrument 
(wavelength bandpass ranging from 5 to 38 A) can resolve strong 
X-ray lines and measure elemental abundances accurately in the 
compact central regions of nearby bright clusters of galaxies. 
We examined the RGS observations of the cool core cluster 
A3 112 (see Table 1 ) using version 1 1.0 of the Science Analysis 


We find that 30.3% ± 5.4% to 37.1% ± 7.1% of the total SN 
explosions are SNe la for models that use W7, CDD, and WDD 
yields. Application of M10 CDDT and CDEF SN la models 
implies a significantly higher fraction of SN la explosions 
(~57% ± 14.5%). The reason for this difference is that M10 


Table 2 

The Best-fit Parameters of the snapec Model Obtained Using T95, 199, and M10 SN Yields and XMM-Newton RGS Spectra 


SNe la 
Model 

SNe cc 
Model 

kT e 

(keV) 

AfSNe 

(xlO 9 ) 

R 

R % 

(%) 

X 2 of Source 
Spectra 
(2608 dof) 

X 2 of Bkgd 
Spectra 
(1490 dof) 

T95 (W7) 

T95 

3.69 ± 0.23 

1.06 ± 0.33 

0.52 ± 0.09 

34.3 ± 6.4 

2817.8 

2026.8 

199 (W7) 

199 

3.70 ± 0.23 

1.06 ± 0.34 

0.52 ± 0.09 

34.2 ±6.3 

2818.7 

2026.4 

199 (W70) 

199 

3.74 ± 0.23 

1.07 ±0.33 

0.50 ± 0.09 

33.5 ± 6.2 

2822.5 

2025.5 

199 (WDD1) 

199 

3.81 ±0.23 

1.15 ± 0.35 

0.57 ±0.10 

36.2 ±6.8 

2835.8 

2022.9 

199 (WDD2) 

199 

3.80 ± 0.22 

1.09 ±0.34 

0.48 ± 0.08 

32.3 ± 5.9 

2831.6 

2023.6 

199 (WDD3) 

199 

3.80 ±0.21 

1.06 ± 0.33 

0.44 ± 0.08 

30.3 ± 5.4 

2829.5 

2024.1 

199 (CDD1) 

199 

3.81 ±0.22 

1.17 ±0.35 

0.59 ±0.12 

37.1 ± 7.1 

2836.6 

2022.8 

199 (CDD2) 

199 

3.81 ±0.22 

1.08 ±0.34 

0.45 ± 0.08 

31.2 ±5.7 

2831.5 

2023.7 

M10 (W7) 

199 

3.70 ± 0.23 

1.06 ± 0.33 

0.51 ±0.09 

33.9 ±6.4 

2819.2 

2826.2 

M10 (CDEF) 

199 

3.66 ±0.30 

1.28 ±0.43 

1.23 ±0.27 

56.3 ± 13.7 

2809.5 

2026.9 

M10 (CDDT) 

199 

3.33 ±0.23 

1.17 ±0.41 

1.34 ±0.29 

57.2 ± 14.5 

2806.4 

2026.8 

M10 (ODDT) 

199 

3.71 ±0.23 

1.12 ±0.34 

0.61 ±0.11 

37.8 ± 7.2 

2823.6 

2025.1 
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Figure 1. First- and second-order RGS sources and background spectral fits using the snapec model. The fit was performed using Iwamoto et al. ( 1 999) SN la yields 
obtained from W7 models and SN cc yields obtained from the weighted average over a Salpeter IMF. The best-fit temperature, total number of supernovae, and SN la 
percentage of the total SNe are 3.69 ± 0.23 keV, (1.24 ± 0.33) x 10 9 , and 34.3% ± 6.4%, respectively (Table 2). The total y 2 is 2818.7 for 2615 degrees of freedom 
for the best fit to the source spectra, and 2026.4 for 1498 degrees of freedom for the background fit. (Table 2). 


CDDT, CDEF SN la models have lower Fe abundance compared 
to the W7, W70, CDD, WDD, and ODDT models. Since Fe is 
mainly produced by SN la explosions, a lower abundance of Fe 
in SN la models in M10 CDDT, CDEF yields results in a higher 
number of SN la explosions and thus higher SN la fraction. The 
fraction of SN la explosions inferred using the 199 CDD and 
WDD models reported in Table 2 are consistent with those in 
de Plaa et al. (2007) based on elements accessible to the EPIC 
detectors: Si, S, Ar, Ca, Fe, and Ni. However, de Plaa et al. 
(2007) find a lower SN la fraction for W7 SN la models. 

The SN models we use in this work produce equally good fits 
to RGS data (see Table 2). Therefore, global fits to the currently 
available high-resolution X-ray observations of clusters of 
galaxies do not allow one to distinguish between different types 
of SN models in an individual cluster. The RGS first- and 
second-order spectral fits obtained using the 199 SN la (W7) 
yields and Salpeter-IMF-averaged 199 SN cc yields are shown 
in Figure 1. We also show a close-up view of the fit in the RGS 
Fe-F bandpass obtained using the T95 SN la (W7) and Salpeter- 
IMF-averaged T95 SN cc yields in Figure 2. Figure 2 shows that 
the T95 SN yields accurately reproduce the RGS spectrum in 
the Fe-F bandpass. 

3.2. XMM -Newton EPIC Observations of A3 112 

We examined the EPIC observations of the cool core cluster 
A31 12 (see Table 1 ) using version 11.0 of the SAS software. The 
EPIC data processing and background modeling were carried 
out with the XMM-Newton Extended Source Analysis Software 
(XMM-ESAS) and methods (Kuntz & Snowden 2008; Snowden 
et al. 2008) described in Paper I. Since the EPIC spectra can 
constrain Si, S, Fe, and O abundances only in the innermost 
30" region (see Paper I), we examined only the central 30" 
spectra in this work. We fit the cluster emission using the 
absorbed snapec model allowing the absorption column density, 
(V SNe , R. redshift, kT e , and normalization parameters freely 
vary. The best-fit parameters of the instrumental and cosmic 
X-ray background model components were fixed to the best-fit 



Figure 2. Close-up view of the Fe-L bandpass. The fit was obtained by applying 
the T95 SN la (W7, slow deflagration) and SN cc (weighted over 10-50 M© 
using the Salpeter IMF function) yields to the XMM-Newton RGS spectra. The 
best-fit model parameters and the goodness of fit are shown in Table 2. 


values obtained from the overall fits and reported in Paper I. 
A conversion factor of 0.41 was again applied to the SNe per 
10 12 M 0 derived from EPIC spectral fits to obtain the /V SXc 
parameter. The best-fit model parameters obtained from the 
snapec fits to EPIC spectra for various SN nucleosynthesis 
models are shown in Table 3. These final yf values, reported 
in Table 3, were obtained from the best fit determined by 
C-statistics. 

The temperature, SN ratio ( R ), and the total number of SN 
measurements obtained from EPIC spectra are consistent with 
RGS results (Table 2) at the lcr level. The lcr difference is due to 
different temperature measurements reported by EPIC and RGS 
observations. Most of the SN nucleosynthesis models we use in 
this work produce equally good fits to EPIC data as in the case 
of RGS analysis. High signal-to-noise CCD resolution EPIC 
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Table 3 

The Best-fit Parameters of the snapec Model Obtained Using T95, 199, and M10 SN Yields and XMM-Newton EPIC Spectra 


SNe la 

SNe cc 

kT e 

AfSNe 

R 

R % 

x 2 

Model 

Model 

(keV) 

(xlO 9 ) 


(%) 

(1683 dof) 

T95 (W7) 

T95 

3.49 ± 0.02 

1.06 ±0.02 

0.40 ± 0.01 

28.6 ± 1.0 

2737.7 

199 (W7) 

199 

3.50 ± 0.02 

1.06 ±0.03 

0.40 ± 0.01 

28.6 ± 1.0 

2737.8 

199 (W70) 

199 

3.50 ± 0.02 

1.06 ±0.02 

0.39 ± 0.01 

28.1 ± 1.0 

2723.9 

199 (WDD1) 

199 

3.47 ± 0.02 

0.82 ± 0.02 

0.76 ± 0.03 

43.2 ± 2.9 

2893.6 

199 (WDD2) 

199 

3.49 ± 0.02 

0.89 ± 0.02 

0.51 ± 0.01 

33.8 ± 1.0 

2783.9 

199 (WDD3) 

199 

3.47 ± 0.02 

0.91 ±0.11 

0.42 ± 0.02 

29.6 ± 1.9 

2749.8 

199 (CDD1) 

199 

3.41 ± 0.01 

0.82 ±0.11 

0.81 ± 0.03 

44.8 ± 2.9 

2927.6 

199 (CDD2) 

199 

3.50 ± 0.02 

0.87 ±0.11 

0.49 ± 0.02 

32.9 ± 1.9 

2767.1 

M10 (W7) 

199 

3.51 ± 0.02 

1.06 ±0.01 

0.40 ± 0.01 

28.6 ± 1.0 

2744.1 

M10 (CDEF) 

199 

3.54 ± 0.02 

1.36 ±0.11 

0.77 ± 0.03 

43.5 ± 2.9 

2775.8 

M10 (CDDT) 

199 

3.39 ± 0.02 

1.11 ±0.10 

1.28 ±0.07 

56.1 ±6.5 

3627.5 

M10 (ODDT) 

199 

3.49 ± 0.02 

0.83 ±0.13 

0.82 ± 0.04 

45.1 ± 3.8 

2842.3 



Figure 3. snapec fit to the XMM-Newton MOS1 (red), MOS2 (blue), and PN 
(magenta) spectra from the innermost 30" region of A3 112 obtained from the 
observation 0603050101. The background model parameters were fixed to the 
best-fit values reported by Bulbul et al. (2012). The fit was performed using 
Iwamoto et al. (1999) W70 SN la and Iwamoto et al. (1999) SN cc yields 
obtained from the weighted average yields over a Salpeter IMF nucleosynthesis 
yields. The best-fit temperature and SN ratio ( R ) and the total number of SN 
(W SNe ) parameters, 3.50 ± 0.02 keV, 0.39 ± 0.01, and (1.06 ± 0.02) x 10 9 , 
respectively, are shown in Table 3. 

(A color version of this figure is available in the online journal.) 

spectra yield stringent constraints on the model parameters (e.g., 
SN ratio ( R ) and the total number of SNe) but cannot distinguish 
between different types of SN nucleosynthesis models. Maeda 
et al. (2010) CDDT SN la models yield a worse fit to the EPIC 
spectra compared to Iwamoto et al. (1999), Tsujimoto et al. 
(1995), or Maeda et al. (2010) W7, CDEF, ODDT SN la models. 

We find that fits to the EPIC spectra find a slightly wider 
range of SN type la fraction, so that 28.1% ± 1.0% to 44.8% ± 
2.9% of the total SN explosions are SNe la for models that use 
W7, CDD, and WDD nucleosynthesis yields. As in the case of 
the RGS analysis, the M10 CDEF, CDDT SN la models imply 
a significantly higher fraction of SN la explosions (43.5% ± 
12.9% to 56.1% ± 6.5%) as a result of lower Fe abundance 
compared to the W7, W70, CDD, WDD models. We also found 
that all W7 models from several simulations (e.g., T95, 199, and 
M10) produce similar goodness of fit, the total number of SNe, 
and SN la fraction. The snapec fits to EPIC MOS1, MOS2, 
and PN spectra performed using the Iwamoto et al. (1999) W70 
SN la and Salpeter-IMF-weighted Iwamoto et al. (1999) SN cc 
metal yields are shown in Figure 3. 


3.3. Astro-H Simulations of A3 112 

The Japan/U.S. Astro-H Observatory is scheduled to launch 
in 2014 and will carry an X-ray calorimeter, the SXS, with a 
high energy resolution of ~4.5 eV (Takahashi et al. 2010). The 
detector will be able to determine line widths and metal abun- 
dances with high precision, revolutionizing our understanding of 
physical processes (e.g., turbulence, the contributions of SN ex- 
plosions to metal enrichment) in the ICM. We simulated 200 ks 
SXS observations of A3 112 using the X-ray events simulator 
software ( simx 1.2.1), based on the temperature and abundance 
measurements obtained from the snapec fits convolved with 199 
SN la (W7) and 199 Salpeter-IMF-weighted SN cc yields, simx 
uses predefined detector response matrix, effective area, and pre- 
dicted background for the SXS instrument to produce an event 
file by convolution with the telescope’s point spread function. 4 
The SXS spectrum is extracted using the FTOOL xselect based 
on the event file produced with the simx software. 

We then fit the simulated SXS spectrum with the snapec 
model convolved with various SN la and SN cc yields for 
variable temperature, /V SNe , R. redshift, and normalization. 
Figures 4 and 5 show the soft X-ray band and a close-up view 
of the Fe-K bandpass for fits to the simulated SXS spectrum of 
A31 12 obtained assuming the same SN yields used to create the 
spectrum. The relative contributions of SN la and SN cc products 
to the spectrum are shown as green and blue lines in Figures 4 
and 5, respectively. The figures illustrate how the significant 
amounts of Fe and Mg produced by SN la explosions, and large 
quantities of O synthesized in SN cc explosions, are manifest 
in the spectrum. The best-fit parameters and goodness of fits 
obtained using the 199, T95, and M10 SN yields are shown in 
Table 4. The features from metals primarily produced by SN 
la explosions such as S, Ar, Ca, and Ni that were inaccessible 
to the EPIC and RGS now become important diagnostics for 
Astro-H data analysis. As a result of the accurate abundance 
measurements of these metals, uncertainties in the /V SNe and 
R parameters are reduced, and they may be determined more 
robustly and with higher precision (see Table 4). 

As expected, good fits are found using SN models that have 
similar elemental yields, e.g., using T95 (W7) or 199 (W70) 
SN la models, to those of the 199 (W7) SN la and SN cc 
yields adopted to produce the simulated Astro-H spectrum (see 
Table 4). On the other hand, the significantly higher x 2 values 
obtained from fits with SN models that have markedly different 


4 http://hea-www.harvard.edu/simx/index.html 
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Energy (keV) 


Figure 4. Simulated Astro-H SXS spectra using the snapec model. The fit was performed using Iwamoto et al. (1999) SN la yields obtained from W7 models and SN 
cc yields obtained from the average yield computed using Salpeter IMF weighting (as in Figure 1). The red color shows the best-fit snapec model with T95 W7 SN la 
and SN cc yields. The contribution of SN la explosions to the metal abundances are shown in green while the contribution of SN cc explosions are shown in blue. 


Table 4 

snapec Parameters for the Best Fit to the Simulated Astro-H Spectra of A31 12 Assuming 199 SN la and SN cc Yields 


SNe la 
Model 

SNe cc 
Model 

kT e 

(keV) 

jySNe 

(xlO 9 ) 

R 

X 2 

(886 dof) 

T95 (W7) 

T95 

3.77 ± 0.03 

1.10 ±0.07 

0.50 ± 0.02 

1124.5 

199 (W7) 

199 

3.77 ± 0.03 

1.09 ±0.06 

0.51 ± 0.02 

1119.9 

199 (W70) 

199 

3.76 ± 0.03 

1.08 ±0.06 

0.49 ± 0.02 

1113.8 

199 (WDD1) 

199 

3.78 ± 0.03 

1.06 ±0.06 

0.60 ± 0.03 

1319.4 

199 (CDD1) 

199 

3.69 ± 0.05 

1.07 ± 0.08 

0.33 ± 0.02 

1348.9 

M10 (CDEF) 

199 

3.76 ± 0.03 

1.46 ±0.08 

0.89 ± 0.05 

1223.9 

M10 (CDDT) 

199 

3.83 ± 0.03 

1.16 ±0.05 

1.54 ±0.12 

1795.0 

M10 (ODDT) 

199 

3.84 ± 0.03 

1.07 ± 0.07 

0.69 ± 0.03 

1292.6 



Energy (keV) 

Figure 5. Close-up view of the Fe-K bandpass of the 100 ks simulated Astro-H 
SXS spectrum, based on the snapec model convolved with 199 SN la and SN 
cc yields. The red color shows the best-fit snapec model for T95 W7 SN la and 
SN cc yields. The contribution of SN la explosions to the metal abundances are 
shown in green while the contribution of SN cc explosions are shown in blue. 

yield patterns, e.g., 199 CDD and 199 WDD models, demonstrate 
that future Astro-H observations of galaxy clusters will allow 
one to distinguish between different SN models and provide 
fundamental constraints on SN nucleosynthesis. 


4. DISCUSSION 

By developing and applying a new spectral analysis method- 
ology to RGS spectra, we measure the total number of SNe that 
have enriched the ICM presently in the inner 52 kpc of A31 12 
and the relative contributions of SNe la and SNe cc. In this sec- 
tion, we compare this with the enrichment expected based on the 
number of SN la and SN cc explosions derived from the amount 
of starlight and well-established empirically measured SN rates 
per star formed. We also investigate prospective local enrich- 
ment sources related to the brightest cluster galaxy (BCG). In 
this way, we can interpret our data analysis in the context of the 
origin of the gas and its metals. 

4.1. Global Cluster Context for the Number of SNe 

Within its virial radius, a sufficiently massive cluster is 
approximately a closed box. For a global efficiency of converting 
gas to stars e S f, the total mass in stars formed is 

4^*form — £sf 44| W y , (10) 

where /V/ hary is the total baryon mass. At the present time the total 
mass in stars, whether contained in individual cluster galaxies 
(including the BCG) or associated with intracluster light (ICL), 
is 

M* = Mbary £sf (1 — R*), (11) 
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1 

S' 

§ 

II 

2? 

( 12 ) 

— ^bary [1 £sf(l ^*)] ? 

(13) 


where R * is the mass return fraction. We can safely neglect the 
distinction between the mass in gas and the mass in the ICM 
for a cluster as rich as A3 112 where the gas mass in galaxies 
is relatively small, and henceforth equate M gas with Micm- The 
star formation efficiency, in terms of the observable M| C m/M*, 

s S f = (1 — + Mjcm/M*)-'. (14) 


The total number 
expressed as 

and 


of SN cc and Type 

n cc = n cc M* form 

A7 Ia = ,/ a M* form , 


explosions can be 

(15) 

(16) 


where r] cc and r] l ' d are, respectively, the specific numbers of SN 
cc and SN la explosions per star formed. Connecting if c and )/ a 
to the snapec model parameters defined in Section 2, 


7 Ia 

R = — , 

r] cc 

(17) 

: 77(1 - R*r'M*, 

(18) 


where rj = rj n + 17 la . We note that these refer to all of the SNe 
while in Section 2; iV SNe and R parameters refer only to SNe 
that enrich the ICM. 

Adopting a “diet Salpeter” IMF that produces relatively fewer 
low-mass stars (Bell & de Jong 2001), Maoz & Mannucci (201 1) 
estimated )/ a to be ~0.002. Adopting the mass return fraction 
and specific SN cc rate for the diet Salpeter IMF: R* ~ 0.35 
(Fardal et al. 2007; O’Rourke et al. 2011), )/ a ~ 0.002, and 
77 11 ~ 0.008 (Maoz & Mannucci 2011; Botticella et al. 2012), 
one predicts R ~ 0.25, which corresponds to Rl'j = 0.2, 
and 77 ~ 0.01. For clusters in the A3112 mass range (e.g., 
Msqo ~ 3.0 x 10 14 Mq, Paper I), recent analysis finds, typically, 
47 */Micm ~ 0. 1 (Balogh et al. 2011; Lin et al. 2012). Therefore, 
from Equation (18), the estimated total number of SN explosions 
per 10 12 M e of ICM is ~1.54 x 10 9 . Multiplying by the 
correction factor of 0.4 1 appropriate to the core of A3 1 1 2, /V SNc 
becomes ~6.3 x 10 8 . Comparing with our results in Table 2, we 
may conclude that the inner 52 kpc of A3 1 12 has been enriched 
by more SNe, with a higher percentage of SNe la, than is true 
globally for a typical cluster of its mass. 


of the class) assuming yields from Kobayashi et al. (2006). This 
enables us to estimate the lock-up corrections, 77 “* and r/ Ia *, to 
the total specific numbers of SNe per star formed. This provides, 
in turn, the corresponding values available to enrich the ICM, 
'I icm an d ! 7 icm : 

/ w Ia * \ 

* = (19) 

= 6.0 X 10 ~ 4 Z Fe *, (20) 

/ N cc * \ 

rr= UrK^ (21) 

= 4.8 x 10“ 3 Z Fe *, (22) 

hicM = >7 Ia - ri la * = 1.4x10 - 3 M- 1 , (23) 

and 

rite m = 7“ - r , | cc * = 3.2 x 10 - 3 M" 1 . (24) 

That is, the metal production from ~60% of SNe cc, and 
~30% of SNe la, must be locked up in stars to enrich them 
to solar Fe abundances and [a/Fe]* ~ 0.25. This results in 
revised predictions of R ~ 0.30 — consistent with our RGS 
results for A3 112, but an even lower expected total number of 
SNe— A SNe ~ 2.9 x 10 8 SNe. 

4.3. Direct SNe Injection 

Sand et al. (2012) recently estimated a specific SN la rate in 
cluster galaxies of rsNia ~ 0.04 SNuM 8 within /Goo (~1 Mpc) 
(see also Maoz et al. 2010) that may be directly injected into the 
ICM by the r < 52 kpc stellar population. The implied level of 
SNIa available for enrichment accumulated over time r is 

N ld = lCT 4 M*(<52kpc) ( — — — ) ( \ , (25) 

v p 4 VO.l SNuM/ Vl0 9 yr7 v 4 

where the stellar mass M* refers to the r < 52 kpc region, and 
TsNia should now be interpreted as an average over r. We use the 
Hernquist approximation (Hernquist 1990) to deVaucouleurs 
profiles for the A3 1 12 BCG and ICL components and estimate 
M* = 3.6 x 10 11 Mq within 52 kpc — i.e., M* = 0.88 Micm- 
A large star-to-gas ratio, M* ~ Micm. and a long timescale 
r > 10 9 yr (with the accompanying increase in r S Nia) would be 
required for this to be a significant source of enrichment. 

4.4. The Cluster Enrichment Puzzle as Realized 
in the Core of A3 112 


4.2. SN Metals Locked Up in Stars 

The galactic mass in clusters is dominated by early-type 
galaxies that form their stars quickly. This results in the well- 
established enhancement in [a/Fe], the abundance ratio of 
a-elements to Fe (expressed as the logarithm with respect 
to solar) — i.e., SNe cc are preferentially locked up in stars 
(Loewenstein 2004). In their investigation of the giant ellip- 
tical galaxy NGC 4472, Loewenstein & Davis (2010) found 
that a ratio of SNe la to total SNe of N u */N SN c * ~ 0.11 
(/V cc *//vSNe* ~ 0.89) and a number of SNe per mass in (present- 
day) stars of 7V SNe /M* ~ 0.0083 resulted in [a/Fe]* ~ 0.25 
and Z Fe * ~ 1 (as observed in this particular galaxy, but typical 


For most of the yield sets we consider, the level of ICM 
enrichment for the central (4.06 ± 0.70) x 10 11 Mq requires 
~1.1 x 10 9 SNe, with SNe la accounting for ~ one-third 
of the total. Although our SN parameters are derived for the 
inner 52 kpc, the high gas mass and temperature imply that 
this gas is primarily intracluster in origin. If we account for 
the fraction of SN products locked up in cluster galaxy stars, 
which is higher for SNe cc based on the measured stellar [a/Fe] 
enhancement, we can explain the proportions of SNe la and SNe 
cc that we infer enriched the ICM. However, the total number 
of SNe we infer is greater than expected by a factor of ~4 for 


5 1 SNuM = Supernova rate per 100 yr per 10 10 A/q in stars. 
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a typical star-to-gas ratio, M.,J M [( - m ~ 0.1. These conclusions 
are altered by direct injection of SN la ejecta in the BCG only 
if accumulated over a large fraction of a Hubble time and are 
exacerbated for the M10 CDEF and CDDT models. Evidently, 
the stellar population responsible for enriching this material is 
characterized by enhanced efficiency in producing SNe, and/or 
the mass in stars responsible for the enrichment, M* > 0.3Micm- 

Larger estimates of ~ 0.2 at the A3 1 12 mass range 

may be found in Gonzalez et al. (2007; see also Langana et al. 
2011). Indeed A3 1 12 was investigated by Gonzalez et al. (2007), 
from which we may derive ~ 0.27. The implication 

is that there is substantial ICL component unaccounted for in 
other studies (see also Bregman et al. 2010) and that the star 
formation efficiency is very high (Equation (14)). 

Even if of intracluster origin, the enrichment of this inner 
cluster gas is likely to reflect the somewhat special conditions 
in the cluster core where high initial overdensity may result 
in enhanced efficiency of galaxy and star formation. Thus, the 
effective value of may very well be larger than is 

representative of the ICM as a whole. At the same time, for rich 
galaxy clusters in general, the global ICM shows enrichment 
beyond what is expected based on the stars we see in galaxies 
today (Portinari et al. 2004; Loewenstein 2006; Maoz et al. 
2010; Bregman et al. 2010). 

Another possible mechanism contributing to the greater than 
expected number of SNe found in our fits is strong early 
enrichment of the lowest-entropy X-ray-emitting gas stripped 
from early galaxies. The heavier metals found in gas stripped at 
large radii will settle into the cluster core due to the central 
gravitational potential, while the stars that produced them 
remain at larger radii. Therefore, a significant fraction of metals 
residing in the cluster core that we observe today may not have 
originated from the local stellar population (Million et al. 2011). 
In this scenario, the apparent total expected number of SNe will 
be enhanced relative to that inferred from A more 

extensive study of a larger sample is needed to fully address this 
issue. 

5. CONCLUSIONS 

In this work we introduce a new approach to determine the 

total number of SN explosions and the relative contribution of 

SNe la to the total enrichment of the ICM using high-resolution 
X-ray observations of clusters of galaxies. The method has been 
implemented as an XSPEC model called snapec. The method 
relies upon the apec model for spectral emission, but uses the 

metal abundances determined from SN yields from the latest 
published SN la and SN cc explosion models normalized to the 

solar abundances (e.g., Anders & Grevesse 1989; Lodders 2003, 

or Asplund et al. 2009). This model returns the total number of 

SN explosions and the relative contribution of the SN types in 
the spectral extraction region that enrich the ICM integrated over 
the cluster life-time, although it does not explicitly constrain the 
SN products locked up in stars in member galaxies. Contrary 
to previous studies in the literature, this method has the ability 
to use all relevant elements, allowing spectra and SN yields 
obtained from published models to determine the best fit to 
their abundances. The main advantage of this approach is that 
it provides a self-consistent single set of uncertainties on the 
total number of SN explosions and relative fraction of SN 
types by directly allowing these parameters to be determined by 
SN yields provided by simulations. This method provides both 
rigorous assessments of uncertainty and improved statistical 


efficiency in estimating the total number of SN explosions and 
relative SN fractions. 

We use XMM-Newton EPIC and RGS observations of A3 1 12 
to validate the method and determine the total number of SNe 
and fraction of SNe la using various published SN yields 
(Tsujimoto et al. 1995; Iwamoto et al. 1999; Maeda et al. 
2010). Since the T95, 199, and M10 SN models produce equally 
good fits to RGS data, current CCD resolution EPIC and high- 
resolution RGS observations of clusters of galaxies do not allow 
one to distinguish between different SN models. The SN la 
fraction reported by EPIC observations are consistent with RGS 
measurements at the 1 a level. We also find that the T95, 199, 
and M10 models produce consistent values of the total number 
of SNe and SN la fractions. 30.3% ± 5.4% to 37.1% ± 7.1% 
of the total SN explosions are SNe la and the total number 
of SN explosions are in the range of (1.06 ± 0.34) x 10 5 * * * 9 to 
(1.28 ± 0.43) x 10 9 * . The best-fit SN la ratios for WDD and 
CDD we obtained from RGS spectra are consistent with de Plaa 
et al. (2007) measurements. However, de Plaa et al. (2007) find 
a lower SN la fraction for W7 SN la models. 

We also applied this method to simulated Astro-H SXS 
observations of A3 112. Astro-H simulations were performed 
using 199 W7 SN la and 199 Salpeter-weighted SN cc models. 
Assuming that 199 W7 models reflects the true nature of the 
ICM enrichment, higher y 2 values obtained from the 199 WDD, 
199 CDD, and M10 CDDT models indicate that we will be 
able to distinguish between different SN models with Astro-H 
observations of clusters of galaxies. 

The estimated mix of SNe la and SNe cc (30% and 70%, 
respectively) expected in the ICM of A3 112 based on well- 
established, empirically measured SN rates per star formed is 
close to what we infer to have enriched the ICM in the inner 
52 kpc of A3 112 based on this new method of high-resolution 
spectral analysis. The inferred level of SN enrichment, ~2.7 x 
10 9 SNe per 1 0 1 2 A77 of ICM, corresponds to ~ 

0. 3-0.4 — compared to ~ 0.1, which has been 

estimated as typical for clusters in the A3 1 12 mass range. 

When applied to RGS observations of a large sample of 
cool core galaxy clusters, this method can provide stringent 
constraints on the contribution of SN explosions to the metal 
enrichment of the ICM and the evolution of the total number 
of SNe la with redshift. This method will also provide an 
independent measure of a key diagnostic of the nature of SN 
la progenitors, the dependence of the evolution of SN la rate to 
the distribution of the delay time between the formation of the 
progenitor system and its explosion as an SN. 

The authors thank Hiroya Yamaguchi for kindly providing 
help on the simx software. We also thank the referee for their 
careful reading of the paper and insightful suggestions. We 
gratefully acknowledge support for this research from NASA 
XMM-Newton grant NNX09AP92G and NASA ROSES-ADP 
grant NNX09AC71G. 
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